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Abstract 

00 . 

q , The particular properties of dynamics are discussed for the dissipatively coupled 

van der Pol oscillators, non-identical in values of parameters controlling the Hopf 
bifurcation. Possibility of a special synchronization regime in an infinitively long 
band between oscillation death and quasiperiodic areas is shown for such system. 
Features of the bifurcation picture are discussed for different values of the control 
parameters and for the case of additional Duffing-type nonlinearity. Analysis of the 
abridged equations is presented. 

Q ! 1. Introduction 

The system of two coupled van der Pol oscillators is the basic model of nonlinear 
dynamics demonstrating the phenomenon of mutual synchronization. It is of interest 
for investigators concerned with synchronization theory and its possible applications (see 
book [1] and refs. below). For example, the system of two coupled van der Pol oscillators 
>• ; was studied by D,G. Aronson at al. [2], D.S. Cohen and J.C. Neu [3, 4] in the case 
of dissipatively coupled oscillators, N. Minorsky [5], R.H. Rand with P.J. Holmes and 
T. Chakraborty [6-8] in the case of weak inertial coupling, R.H. Rand and D.W. Storti 
[9] and Pastor - Diaz at al. [10] in the case of strong inertial coupling. Similar system is 
used also in works of T. Pavlidis [11], M. Poliashenko [12, 13] and D.S. Cohen, J.C. Neu 
[4, 5] for modeling biological and chemical processes. Description of the system by means 
of Adler phase equation was presented in work [7]. However, nowadays interest to this 
problem does not decrease, since new aspects and oscillation effects are detected. In this 
respect we may point out, for example, recent resumptive work [14], where this system 
is studied within the framework of quasiharmonic approximation including possibility of 
combined coupling. However, most of the works on this subject does not go beyond 
the case of identical oscillators' parameters controlling the Andronov - Hopf bifurcation. 
Work [14] is an exception, but the special attention is also devoted there to the case 
of identical subsystems. Moreover, authors have restricted themselves to the assumption 
that controlling parameters are small, also they didn't take into account the anisochronism 
of oscillators. But it is anisochronism that results in the normal form of the Andronov - 
Hopf bifurcation and its consideration is important in respect to the generalization of 
results. 

In this work we shall carry out the investigation of the parameter space structure 
and of possible oscillation regimes in the dissipatively coupled van der Pol and van der 
Pol - Duffing oscillators with non-identical controlling parameters. The initial system of 
differential equations describing the interaction between the oscillators is of the form 

d oc , - o > dx 'j . doc dti . 
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Here Ai and A 2 are parameters characterizing the excess above the threshold of the 
Andronov - Hopf bifurcation in autonomous oscillators, 6 is the frequency mismatch be- 
tween the autonomous second and first oscillators and \x is the coefficient of dissipative 
coupling. Parameter (3 is related to the insertion of nonlinearity of Duffing oscillator type 
into the van der Pol equation. Within the framework of quasiharmonic approximation 
it is responsible for the nonlinearity in the phase equation and, correspondingly, for the 
anisochronism [1]. Strictly speaking, anisochronism is also possible in the system (1) 
for large enough parameters Ai and A 2 values. But nevertheless we shall refer to the 
parameter (3 as to the parameter of phase nonlinearity or anisochronism. 

As we have already noticed, parameters Ai and A 2 are responsible for the Andronov - 
Hopf bifurcations in subsystems. They define the sizes of limit cycles of autonomous 
oscillators in the phase space. If the values of Ai and A 2 are small, the limit cycles are 
small too. In that case one can use quasiharmonic approximation. We shall choose the 
values of Ai and A 2 from the [1,2] interval, i.e. Ai, A 2 G [1,2]. Such choice gives, on the 
one hand, an opportunity to reveal features that are characteristic for the quasiharmonic 
approximation on the qualitative level at least and, on the other hand, to illustrate the 
phenomena intrinsic to the system (1) that appeared to be essentially richer. 

2. The structure of the (frequency mismatch - coupling value) parameters 
plane for non-identical subsystems. The possibility of the broadband syn- 
chronization 

At first, we shall present the results of computer investigation of the system (1) pa- 
rameter space. We shall use the method of dynamic regime chart construction (see, e.g., 
[15, 16]). Within the framework of such a method we shall mark the oscillation period 
of the system of coupled oscillators by means of dark color and gray color hue on the 
parameter plane (frequency mismatch 5 - coupling value //). White color corresponds 
to the chaotic or quasiharmonic motions. Cycle periods were calculated by means of 
the Poincare section method: this is the number of points of intersection of the phase 
trajectory on the attractor and the surface selected as the Poincare section. Only those 
crossings were taken into account that correspond to the trajectories coming to the surface 
from the one side. 

The system under investigation is characterized by four-dimensional phase space 
(x,x,y,y). Therefore three-dimensional hypersurface that is preset by means of some 
additional condition, e.g., zero velocity of the second oscillator y = may serve as the 
Poincare section. In that case the number n of points of intersection of trajectory and 
section was determined. Colors on the charts are chosen in accordance with the period n. 

The chart of dynamic regimes obtained in such a way for identical van der Pol oscilla- 
tors is given in fig. 1 on the (5, //)-plane for Ai = A 2 = 1. Subsequent characteristic areas 
shown on this chart are: 

• main synchronization tongue with the frequency ratio of 1/1; 

• area of quasiperiodic regimes with the embedded system of higher synchronization 
tongues among which the tongue with the rotation number of 1/3 is the most cha- 
racteristic; 

• area of the "oscillation death" effect [1, 2] which corresponds to the stability of the 
equilibrium state point at the origin in the presence of sufficiently strong dissipative 
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We will characterize mutual oscillations of oscillators by means of the rotation number 
w. By analogy with [16, 17] we find numerically the average return time r y for the chosen 
Poincare section y = and the average return time r x for the section of the first oscillator 
defined by x = 0. Then the rotation number is defined as w = T y /r x . The graph of the 
rotation number w dependency on the frequency mismatch 5 is shown in fig. 2 (a) for the 
coupling parameter value /x = 0.5. On this graph one can see characteristic wide "steps" 
that correspond to the rotation number values w — 1, 1/3, 1/5, etc. 

Discrepancy between the controlling parameter values of coupled van der Pol oscilla- 
tors leads to the change of the (frequency mismatch - coupling value) parameter plane 
structure. The chart of dynamic regimes for non-identical subsystems with Ai = 1.25 and 
A2 = 1 is displayed in fig. 3. Non-identity is important but it is still not very large. Phase 
plane portraits computed at several chosen points of the parameter plane are also shown 
in fig. 3 in pairs on (x, x) and (y, y)-planes. 

Comparing the charts of regimes in fig. 1 and fig. 3, one can see that they are essentially 
different in one respect. The boundary between the oscillation death area and the area of 
quasiperiodic regimes in fig. 3 is not a line but a band of finite width in coupling parameter 
A 2 < /i < Ai that stretches infinitely into the area of increased coupling mismatch. 
Existence of synchronization in the system in the presence of arbitrarily large values of 
oscillators eigenfrequency mismatch one can call "broadband synchronization". 

Moving inside this band towards the increase of the second oscillator's eigenfrequency 
one can observe synchronous regimes that differ from each other in "twistedness" of the 
phase trajectory on the phase plane portraits of the second oscillator presented in fig. 3. 
Corresponding chart areas look like "degraded" tops of synchronization tongues (it is 
mainly visible for the characteristic synchronization tongue w = 1/3). 

The graph of the rotation number w dependence upon the frequency mismatch S inside 
the area of the "broadband synchronization" is presented in fig. 2 (b) for Ai = 1.25, A 2 = 1 
and /i = 1.05. One can see considerable alterations in comparison with fig. 2 (a). Areas 
with irrational rotation numbers disappeared and the system of steps with w = 1,1/3 
and 1/5 became essentially more evident. 

The dynamic regimes chart for the system of coupled van der Pol oscillators (1) is 
shown on the top piece of fig. 4 for the case of distinctly different controlling parameters 
of subsystems Ai = 2, A 2 = 1. As an illustration of the second oscillator's behavior the 
examples of its phase plane portraits on the (y, y)-plane, Poincare sections and y{t) - 
realizations are given on fig. 4. We should remind that the Poincare section was chosen 
as the hypersurface y = in four-dimensional phase space (x, x, y, y). Hence the Poincare 
map is automorphism of the three-dimensional space (x, x, y) as is shown on fig. 4. One 
can see fixed points and cycles with periods of 1, 3 and 5 which correspond to the regime 
types on the chart. The graph of the dependence of the corresponding rotation numbers 
upon the frequency mismatch is shown in fig. 2 (c) for Ai = 2, A 2 = 1 and /i = 1.5. One 
can see steps with w — 1, 1/3 and 1/5 which are changing each other successively. 

The dynamic regime chart for the system of coupled van der Pol oscillators is shown 
on fig. 5 for the case of even greater value of the controlling parameter of the first system 
Ai = 3, A 2 = 1. In contrast to fig. 3, oscillators are essentially non-identical. As fig. 5 
shows, besides the /i-directionally broadening of the synchronization band the considerable 
broadening of domains with fractional rotation numbers inside the band takes place during 
the increase in difference of controlling parameters. These areas are filling gradually 
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the most part of the band of broadband synchronization. So the remarkable feature 
of the problem is the preservation of fixed rotation number in essentially broad range 
of frequency which is several times more than the width of corresponding tongues for 
identical subsystems. 

In case of distinctly different Ai and A2 parameters values one can provide an explana- 
tion of broadband synchronization appearance. It is based essentially on the non-identity 
of oscillators. Really, if /x exceeds both Ai and A2, both oscillators are behind the threshold 
of the "oscillation death" effect. But in the range of A2 < [i < Xi only the second oscillator 
appears to be essentially dissipative. And at the same time the first oscillator appears to 
be leading and excites, in fact, the second one. In this respect, different scales along the 
coordinate axes on the phase plane portraits for the first and the second oscillators are 
quite characteristic. 

To discuss the mutual influence of coupled van der Pol oscillators in the described 
above domain of the parameters plane, we consider the system (1) with the "shut off" 
influence of the second oscillator upon the first one 

d , o , dx 
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The chart of dynamic regimes for the system (2) is given on fig. 6 for Ai = 2, A 2 = 1. 
Comparing this chart with the chart for the case of mutual influence of oscillators (see 
fig. 4) one can see that formation of the broadband synchronization strip is the result of 
the influence of the first oscillator that is behind the threshold of the Andronov - Hopf 
bifurcation upon the second oscillator, oscillations of which would damp without such an 
influence. Size and structure of the broadband synchronization area on fig. 4 and fig. 6 are 
similar. The only difference is that without the mutual influence the oscillation death is 
observed at any point of parameter space satisfying X± < fj, and A 2 < fJ>. In case of absence 
of the second oscillator's influence upon the first one this fact has simple explanation: the 
first oscillator appears to be before the threshold of the Andronov - Hopf bifurcation 
taking place when Ai — /j, — and its oscillations are damping. The second oscillator is 
also below the threshold of the Andronov - Hopf bifurcation, because \i > X\ > A 2 , and 
its oscillations are also damping without the excitation of the first oscillator. 

However, our consideration is valid only when Ai and A 2 values are essentially different. 
If these parameters values are close, this argumentation can't be used. 

3. The analysis of the broadband synchronization by means of the abridged 
equations 

It is interesting to find out, which elements of outlined picture one can reveal within 
the framework of quasiharmonic approximation. We define for it 

x = I ( ae « + a *e- u ) , y = - {be u + b* e - lt ) , (3) 

and use the traditional additional assumption of 

\{ae lt + a*e- u ) = 0, hbe u + 6*e _i *) = 0. (4) 
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After the averaging we obtain equations for the complex amplitudes 
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After the transformation of variables: r = t/2, z = a/2, uj = b/2 the equations (5) yield 

(6) 



dz 

— = Xiz — \z\ 2 z + 3if3\z\ 2 z + fj,(u — z), 
dt 

—— = \ 2 uj — \uj\ 2 uj + 3if3\uj\ 2 u + fj,(z — uj) + ibuj. 
dt 

Setting z(t) = R(t) exp(iipi) and uj{t) = r(t) exp(i</? 2 ) we obtain equations for the 
amplitudes R, r and their phases ipi and <p 2 

dR 

— = R(Xi - fi) - R 3 + fir cos(y? 2 - ipi), 
dr 

dv 

— = r (A 2 - /i) - r 3 + uRcoafa - </? 2 ), 
^ = 3/5i? 2 + ^sin(v9 2 -^i), 

— — = Apr H /j, sm((/?i - </? 2 J + 0- 
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There is only the phase difference in the equations for amplitudes. Therefore, one can 
subtract the second equation for the phase from the first one and obtain the following 
abridged equations containing the relative phase of oscillators ip = <fi — y? 2 : 

dR 

— — = R(\\ — /i) — R 3 + ixr cos ^ 
dr 

dv 

— = r (X 2 -fJ.) -r 3 + fiR cos ip, (8) 
dr 

^ = -5 + 3(3(R 2 - r 2 ) - p(L + -) sin^. 
dr R r 

We have obtained the system of equations similar to [1] but it is taking the possible 
non-identity of controlling parameters into account. 

Evaluate the lower boundary of the broadband synchronization area analytically using 
the abridged equations (8). Consider the stationary system of (8): 

= R(Xi — n) — R 3 + fxr cos -0, 

= r(A 2 - fx) - r 3 + /iRcosi/j, ^ 

S = 3P(R 2 -r 2 )-fi(^ + -)smij. 

R r 

We suppose that the summands fiR cos ip and /ir cos ip are small on the boundary of 
the synchronization area. Then the evaluations of the steady amplitudes of limit cycles 

are 

R w a/Ai - /i, r ~ a/A 2 - /i. (10) 
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Since radical expressions should be non-negative it is necessary to satisfy the conditions 
(Ai > fi) and (A 2 > fi). Substituting (10) into the third equation of (9), we obtain 

6 = 3/?((v^7^) 2 - (v/AT^) 2 ) - li (^Eji + V^H^ sin ^. (ll) 

\VAi-A 4 V^-I^J 

This equation is a sort of the well-known stationary Adler equation [1]. Its stability 
threshold is defined by the condition of sin ^ = 1. But then we have cosip = and our 
assumption of the neglected terms smallness in the first and the second equations (9) is 
justified. 

Consider the case of absence of the phase nonlinearity. Then (11) yields 

5 = ±J ^El + ^A. (12) 



5(/x)-function given by the relation (12) defines the boundary of the synchronization area 
which ought to move to the lower boundary of the infinite synchronization band according 
to fig. 3 and fig. 4. Denote this branch as 6 +([/,) and consider its asymptotic properties 
in case of Ai > A 2 . If the coupling parameter tends to the value of controlling param- 
eter of the second oscillator, then (A 2 — /x) — > 0. It means that \J (A 2 — A t )/(Ai — /x) <C 
a/ (Ai — /u)/(A 2 — /x) and <5(/x)-function tends to the infinity, i.e. lim M _^ 2 <5+(a0 = oo. So 
it must be \i — > A 2 when 5 — > oo for the graph of inverse function \i = /i(5). It means 
that the boundary of the synchronization area extends arbitrarily far to the area of large 
frequency mismatches. 

5 (fi) -function given by the relation (12) defines the lower boundary of the infinitely 
long synchronization band. The upper boundary of this band is the boundary of the 
oscillation death area when the trivial rest state of point of origin becomes stable. Use 
the abridged equations obtained before for the system of coupled van der Pol - Duffing 
oscillators (6) for the analytical estimation of this boundary. After the linearization of 
the system (6) about a point of origin we obtain 

dz . , . 

— = \ t z + n(u - z), 

%> > r , , <13) 

— — = A 2 uj + U(z — uj) + toco, 
dr 

For the investigation of the stability of zero rest state of the obtained linear system set 

where rj, e are real numbers. It is necessary to satisfy the condition 77 < for the realization 
of the oscillation death regime. Correspondingly, there is r\ = on the boundary of the 
oscillation death area. Then 

z ~ e l£ \ uj ~ e l£t . (15) 
Substitution of (15) in the system (13) yields the relation 



(is — Ai + (J,)(ie — A 2 + ji — i8) = fi 2 <^ 
s 2 + e(-<J+i(Ai + A 2 - 2fj)) + ( f i(X 1 + A 2 ) - AiA 2 + i8(n - Ai)) 



Denote 



ot x = -5, 

a y = \ 1 + \ 2 - 2fi, ^ 
Px = /i(Ai + A 2 ) - AiA 2 , 
(3 y = 6(fj, - Ai). 

Then 

e 2 + e(a x + ia y ) + + i(3 y ) = 0, (18) 

where a x , a y , (3 X , (3 y are real numbers. Therefore one can divide the equation (18) into the 
real and imaginary parts: e 2 + ea x + f3 x = and ea y + f3 y = 0. It follows that 

e = -4 4-^ + ^ = 0. (19) 

a w a tf a w 

Substituting (17) in (19), we obtain after some transformations the relation that defines 
the boundary of the oscillation death area and the upper boundary of the synchronization 
band, correspondingly 

A 2 _ (/j(Ai + A 2 )-A 1 A 2 )(A 1 + A 2 -2^) 2 

(Ai-/i)(A 2 -/i) • [2U) 

In special case when Ai = A 2 = A the relation (20) yields 

5 2 = 4A(2/i - A), (21) 
and the oscillation death area is defined by the certain [1, 2] two-sided inequality 

a< "4( a+ £)' (22) 

Investigation of the behavior of <5(/z)-function given by the relation (20) yields its 
asymptotic features 

lim 5(fj) = oo, (23) 

lim 8(n) = 2 y / f j,(X 1 + X 2 ). (24) 

The graphs of 5(/z)-functions given by the relations (12) and (20) are shown in fig. 7. 
Inside the grey area the rest state of the Adler equation exists and, correspondingly, the 
synchronization of oscillators is possible. We see that there is a synchronization band, 
which is infinitely long in eigenfrequency mismatch, in the coupling parameter domain of 
A 2 < /i < Ai. Consequently, the analysis that was made above allows to show the presence 
of the infinite on frequency mismatch synchronization band. It should be noticed that 
this analysis does not allow to diagnosticate the multiple synchronization regimes and, 
therefore, it is a question of the synchronization of 1/1-type. 
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4. The broadband synchronization in system with the phase nonlinearity 

We consider now the case when there is the phase nonlinearity presented in the sys- 
tem. The question of interest is how the phase nonlinearity influences the broadband 
synchronization domain and the overall structure of the parameter plane (eigenfrequency 
mismatch - coupling strength) . One can determine the lower boundary of the synchro- 
nization band from the equation (11) using the fact that this equation is the stationary 
Adler equation. This boundary is given by the relation 

5 = ±» (^pS + ^SEB + miV^f - (a/A^) 2 )- (25) 

The graph of the 5 (//) -function given by the relation (25) is shown in fig. 7 by the 
dashed line. At the same time, the upper boundary of the synchronization band shown in 
fig. 7 is determined, as before, by the relation (20) since it corresponds to the linearized 
system and is valid independently of the presence or absence of the phase nonlinearity. 
One can see that presence of the phase nonlinearity leads to the displacement of the 
synchronization tongue along the axis of eigenfrequency mismatch and does not influence 
on the broadband synchronization area in the case of large mismatches. 

Dynamic regime charts obtained via numerical simulations of the system of coupled 
van der Pol - Duffing oscillators are shown in fig. 8 for different values of Ai and A2 
in presence of the phase nonlinearity in system. One can see that the synchronization 
tongue shifts to the right along the axes of eigenfrequency mismatch with the increase 
of the parameter of nonlinearity f3, just as the analysis of the Adler equation predicts. 
One can see though that the phase nonlinearity influences strongly on the areas with 
fractional rotation numbers. Inside the synchronization band they decrease slightly in 
size, change their own shape and their frequency threshold shifts strongly to the area of 
large mismatches (to the right). 

Phase plane portraits for different points of the parameter plane are given in fig. 8 (b) 
in the (x, x)-plane (left one) and in the (y, y)-plane (right one). Moving along the syn- 
chronization band, one can see the transition from the synchronization regime of 1/1-type 
to the regime of 1/3-type. One can note one more feature of the phase plane portraits in 
the area of 1/3 in fig. 8 in comparison with the fig. 3. It is the fact that "loops" of the 
second oscillator's attractor on the phase plane portrait may enclose a coordinate origin. 



5. System with non-identity in value of nonlinear dissipation 

In described picture of interaction between oscillators the first oscillator dominates 
more over the second one with increase of the frequency mismatch in the broadband 
synchronization area, when X 2 < /i < Xi- It becomes apparent in noticeable downsizing of 
the second oscillator's limit cycle in relation to the size of the first oscillator's limit cycle. 
(Pay attention to the scales in coordinate axis on the phase plane portraits in fig. 3). But 
it is possible to obtain synchronization regimes when oscillators are equivalent also in case 
of large frequency mismatch. Consider for that the following system of equations: 
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Here additional parameter 7 characterizes nonlinear dissipation of the second oscillator. 
This parameter together with the parameter A defines the size of the limit cycle in au- 
tonomous behavior. Using quasiharmonic approximation one can estimate corresponding 
sizes of the limit cycles: for the first oscillator it is a/Ai and for the second one it is 
a/ A2/7. So using small values of parameter 7 one can greatly enlarge the size of the 
second oscillator's limit cycle. 

Dynamic regime chart for the system (26) for Ai = 2, A2 = 1,7 = 0.01,(3 = and 
several typical phase plane portraits are shown in fig. 9. One can see that the broadband 
synchronization remains in that case. Moreover, the main synchronization tongue appears 
to have an extension within which the second oscillator dominates over the first one (see 
fig. 9 (a)). In areas of multiple synchronization it may be both the situation of equivalent 
oscillators and the situation of changing of "leading" oscillator, when the second oscillator 
dominates over the second one (see fig. 9 (b)). "Equality" or oscillator's domination may 
be determined from the shape of the phase plane portrait and the size of the attractor. 

6. Conclusions 

Dissipatively coupled non-identical in controlling parameters van der Pol oscillators 
demonstrate interesting synchronization features. On the parameter plane (frequency 
mismatch - coupling value) the boundary between the oscillation death area and the area 
of quasiperiodic regimes turns from the line into the band of finite width which stretches 
along the mismatch parameter axis into the area of large mismatch values. In the case 
of small non-identity the synchronization with the rotation number of 1/1, generally, 
corresponds to this band. However, there are also regimes of 1/3, 1/5 types, etc. and their 
domains of stability on the parameter plane look like "degraded" tops of corresponding 
synchronization tongues. Areas corresponding to these regimes expand essentially with 
the increase of non-identity, filling, practically, the whole band. The oscillation regime 
when one oscillator dominates distinctly over another oscillator is realized in that case. 
Areas with fractional rotation numbers decrease in size and their frequency threshold 
moves towards the large mismatch values in the case of presence of the phase nonlinearity 
of Duffing oscillator type. The size of the "broadband" synchronization area increases 
with insertion of a non-identity in value of nonlinear dissipation due to the formation of 
"outgrowth" from the main synchronization tongue. And at the same time regimes, for 
which oscillators are equivalent and the changing of "leading" oscillator may take place, 
appear. 

Acknowledgements 

This work was supported by RFBR (Project No. 06-02-16773) and foundation of 
nonprofit programs "Dynasty". 

References 

[1] A. Pikovsky, M. Rosenblum, J. Kurths [2001] Synchronization (Cambridge), p. 411 

[2] D. G. Aronson,G. B. Ermentrout, N. Kopell [1990] "Amplitude response of coupled 
oscillators", Physica D41, 403-449. 



n 



[3] D. S. Cohen, J. C. Neu [1979] "Interacting oscillatory chemical reactors", in Bifurca- 
tion Theory and Applications in the Scientific Disciplines, eds. Gurel, O. and Rossler, 
O. E. (Ann. N.Y. Acad. Sci. 316) pp. 332-337. 

[4] J. C. Neu [1979] "Coupled chemical oscillators", SI AM J. Appl. Math. 37(2), 307-315. 

[5] N. Minorsky [1962] Nonlinear oscillators (Van Nostrand). 

[6] R. H. Rand, P. J. Holmes [1980] "Bifurcation of periodic motions in two weakly 
coupled van der Pol oscillators", Int. J. Non-Linear Mechanics 15, 387-399. 

[7] T. Chakraborty, R. H. Rand [1988] "The transition from phase locking to drift in a 
system of two weakly coupled van der Pol oscillators", Int. J. Non-Linear Mechanics 
23(5/6), 369-376. 

[8] T. Chakraborty [1986] "Bifurcation analysis of two weakly coupled van der Pol os- 
cillators", Doctoral thesis (Cornell University). 

[9] D. W. Storti, R. H. Rand [1982] "Dynamics of two strongly coupled van der Pol 
oscillators", Int. J. Non-Linear Mechanics 17(3), 143-152. 

[10] I. Pastor-Diaz, A. Lopez-Fraguas [1995] "Dynamics of two coupled van der Pol oscil- 
lators", Phys. Rev. E52, p. 1480. 

[11] T. Pavlidis [1973] "Biological oscillators: the mathematical analysis", Academic 
Press. 

[12] M. Poliashenko, S. R. McKay, C. W. Smith [1991] "Chaos and nonisochronism in 
weakly coupled nonlinear oscillators", Phys. Rev. A44, p. 3452. 

[13] M. Poliashenko, S. R. McKay, C. W. Smith [1991] "Hysteresis of synchronous - 
asynchronous regimes in a system of two coupled oscillators", Phys. Rev. A43, p. 
5638. 

[14] M. V. Ivanchenko, G. V. Osipov, V. D. Shalfeev, J. Kurths [2004] "Synchronization 
of two non-scalar-coupled limit-cycle oscillators", Physica D189, 8-30. 

[15] A. P. Kuznetsov, S. P. Kuznetsov, I. R. Sataev [1997] "A variety of period-doubling 
universality classes in multi-parameter analysis of transition to chaos", Physica 
D109, 91-112. 

[16] A. P. Kuznetsov, E. Mosekilde, L. V. Turukina [2006] "Synchronization in systems 
with bimodal dynamics", Physica A121(2), 280-292. 

[17] D. E. Postnov, A. V. Shishkin, O. V. Sosnovtseva, E. Mosekilde [2005] "Two-mode 
chaos and its synchronization properties", Phys. Rev. E72, pp. 056208(5). 



1 n 




oscillation 
death area 



/ 



Fig. 1: Dynamic regime chart for the system (1) for Ai = A2 = 1, j3 = 0. 
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Fig. 2: Dependence of the rotation number w upon the frequency mismatch 5 for 
(a) Ai = 1.25, A 2 = 1, n = 0.5;(b) Ai = 1.25, A 2 = 1, n = 1.05; (c) A x = 2, A 2 = 1, 
^ = 1-5. 
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Fig. 3: Dynamic regime chart for the system (1) for Ai = 1.25, A 2 = 1, (3 — and phase 
plane portraits in characteristic areas of the parameter plane. 
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Fig. 4: Dynamic regime chart for the system (1) for Ai = 2, A2 = 1, f3 — and phase 
plane portraits of the second oscillator, its Poincare sections and realizations in typical 
areas of the parameter plane. 
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Fig. 5: Dynamic regime chart for the system (1) for Ai = 3, A2 = 1, (3 — 0. 
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Fig. 6: Dynamic regime chart for the system (2) with the "shut off" influence of the second 
oscillator upon the first one for for Ai = 2, A 2 = 1. 



1 c 




Fig. 7: Synchronization area of 1/1-type defined analytically within the framework of the 
phase approximation for Ai = 2, A2 = 1, (3 — shown by the grey color. Boundary of the 
synchronization area in presence of the phase nonlinearity f3 — 1 is shown by the dashed 
line. 
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Fig. 8: Dynamic regime charts for the system (1) for Ai = 2, A2 = 1 and (a) (5 = 0.5; 
(b) (5 — 1 with phase plane portraits at characteristic points. 
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Fig. 9: Dynamic regime charts for the system (26) for Ai = 2, A2 = 1, 7 = 0.01, (3 — and 
typical phase plane portraits of the first (on top) and the second (from below) oscillators. 
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